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In this paper we present an overview of radio frequency (RF) spectroscopy in the atomic Fermi superfluids. 
An ultimate goal is to suggest new directions in the cold gas research agenda from the condensed matter perspec- 
tive. Our focus is on the experimental and theoretical literature of cold gases and photoemission spectroscopy 
of the cuprates particularly as it pertains to areas of overlap. In addition to a comparison with the cuprates, this 
paper contains a systematic overview of the theory of RF spectroscopy, both momentum integrated and momen- 
tum resolved. It should be noted that the integrated and momentum resolved forms of photoemission are equally 
important in the high Tc cuprate literature. We discuss the effects of traps, population imbalance, final state in- 
teractions over the entire range of temperatures and compare theory and experiment, most notably in the context 
of recent tomographic scans in population imbalanced gases. We show that this broad range of phenomena can 
be accomodated within the BCS-Leggett description of BCS-BEC crossover and that this scheme also captures 
some of the central observations in photoemission experiments in the cuprates. In this last context, we note that 
the key themes which have emerged in cuprate photoemission studies involve characterization of the fermionic 
self energy, of the pseudogap and of the effects of superconducting coherence (in passing from above to below 
the superfluid transition temperature, Tc). These issues have a counterpart in the cold Fermi gases and it would 
be most useful in future to use these atomic systems to address these and the more sweeping question of how to 
describe that anomalous superfluid phase which forms in the presence of a normal state excitation gap. 

PACS numbers: 



I. INTRODUCTION AND MOTIVATION 

There is considerable excitement surrounding the discov- 
ery fi, 2, 3, 4, 5, 6, 7, 8, 9, lOJ of superfluidity in the ultracold 
Fermi gases. What is novel about these new superfluids is 
that one can tune the attractive interaction from weak (as in 
the BCS limit) to strong as in the Bose Einstein condensation 
(EEC) regime. These experiments will continue to impact 
condensed matter physics by providing, at the least, a new 
class of "materials" which elucidate a very powerful general- 
ization of BCS theory. A number of people HI [H [H H 
have also argued that this BCS-BEC crossover might be rele- 
vant to the cuprate superconductors. Because of their anoma- 
lously short coherence length it is claimed [ 15] that these ma- 
terials are "mid-way between BCS and BEC". That is, the 
attractive interaction driving the superconducting pairing may 
be stronger than that in conventional superconductors. In this 
way the tuneability of the interaction strength in the Fermi 
gases provides an ideal model system with which to study the 
physics of the short coherence length cuprates and the role 
of strong attraction (generally associated with high transition 
temperatures). From a very different perspective, it has also 
been argued that in future optical lattice experiments lfT6ll in- 
volving the atomic Fermi gases, one will be able to simulate 
repulsive Hubbard models and thereby investigate the "Mott 
physics" aspects [17.1 of high Tc superconductivity. 

While condensed matter physicists have a wealth of well- 
developed techniques for characterizing electronic supercon- 
ductors, the tools currently available to the atomic physicists 



who study the Fermi gases are far more limited. Moreover, it 
is not at all straightforward to determine something as com- 
monplace as the temperature in the gas, although some im- 
pressive progress jsl ll Si ll 9tl has been made along these lines. 

This paper is devoted to addressing one of the most power- 
ful techiques currently being applied to the Fermi gases: radio 
frequency (RF) spectroscopy. We will show how this tech- 
nique is similar to that of photoemission in condensed mat- 
ter physics and exploit the analogy, already discussed in the 
literature |20], between momentum resolved RF and angle 
resolved photoemission spectroscopy (ARPES). As a back- 
ground for both communities, we review some of the experi- 
mental and theoretical literature on RF spectroscopy (of cold 
gases) and photoemission spectroscopy (of the cuprates). We 
argue that there are a number of issues which have been cen- 
tral to high temperature superconductivity which would be 
useful to address more systematically in the ultracold Fermi 
gase s. Perhaps the most notable example of commonality 
il ll [T2I1 in this regard is the ubiquitous pseudogap phase 
which is at the core of current studies in the high Tc supercon- 
ductors and has emerged as important in the ultracold Fermi 
gases. 

We begin by focusing on the overlap of the experimental 
concepts behind photoemission experiments ll2ll I22I1 in the 
cuprates and RF spectroscopy in the atomic Fermi gases. We 
will see that both experiments reflect the behavior of the all 
important fermionic spectral function A(k, uj) which charac- 
terizes completely the single fermion or one-particle prop- 
erties of a given many body system. In simplistic terms, 
the driving force motivating the photoemission studies in the 
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cuprates is to acquire an understanding of the "mechanisms" 
and nature of superconductivity. There has been a recent em- 
phasis on high temperatures near T*, where the pseudogap 
turns on and on the region from above to below the superfluid 
transition temperature, Tc- By contrast in the ultracold gases, 
the RF spectra have been used to characterize the pairing gap 
A- much like tunneling is used in conventional superconduc- 
tors. There has been a recent emphasis on very low tempera- 
tures T « Tc and in particular in quantifying the size of A 
at T = 0. 

Some of the key issues which have emerged in photoemis- 
sion studies of the cuprates involve (i) a characterization of 
the self energy contained in the spectral function. Different 
empirical models [23] have been deduced which, it is argued, 
might ultimately hold the clue as to the nature of the mediat- 
ing boson, (ii) Also important is the origin of the all important 
pseudogap. There is a debate ] 1 1, 17] about whether this gap 
is a signature of a hidden order parameter or whether it reflects 
the incipient pairing which ultimately leads to the condensed 
phase at lower T. (iii) It is viewed as extremely important to 
arrive at an understanding of how superconducting coherence 
manifests itself in this spectroscopic experiments as one goes 
from the normal to the ordered phase. This is a complicated 
question, given the presence of a normal state (pseudo)gap. 
Finally, other issues of interest are the nature of the order pa- 
rameter and pseudogap symmetry (which have been shown to 
be consistent with d-wave). 

In the cold gases an underlying goal has been to test differ- 
ent theories of BCS-BEC crossover, particularly establishing 
the most suitable ground state and its quantitative implications 
such as the pair size ].24,] . The parameters which quantify 
the nature of the scale-free or "unitary" gas have also been 
addressed. Of additional interest are studies on how popula- 
tion imbalance ]25, 26, 27, 28] can co-exist with superfluidity. 
Here new phases associated with, for example, the exotic ll29ll 
Larkin- Ovchinnikov, Fulde-Ferrell (LOFF) form of pairing 
have been contemplated. Even more topical is the behavior in 
the limit of extreme imbalance 0271 12811 . 



One can see that, despite the similarities in these two spec- 
troscopic techniques, the research agenda in the two commu- 
nities is rather different. In the high temperature supercon- 
ductors, the focus has been around the temperature regime 
near Tc- Furthermore, quantitative issues are viewed as of 
considerably less importance than aiTiving at a qualitative un- 
derstanding, which is still very incomplete. By contrast in the 
ultracold Fermi gases the focus has been on temperatures as- 
sociated with the ground state and on arriving at a more com- 
plete quantitative characterization. 

This brings us to a major goal of the present paper which 
is to suggest new directions in the cold gas research agenda 
from the condensed matter perspective. In particular we wish 
to highlight differences and similarities in the cold gases with 
the analogous cuprate studies. A general theme, which takes 
a cue from the copper oxide superconductors, is to focus on a 
characterization of (i) the fermionic self energy, (ii) the pseu- 
dogap phase and (iii) how superfluid coherence is established 
and manifested (in these spectroscopies) at and below Tc- 



A. Comparing and Contrasting RF Witli Photoemission 

Photoemission and Angle Resolved Photoemission Spec- 
troscopy (ARPES) have been remarkable tools for charac- 
terizing the cuprate superconductors Here one in- 
vokes the "sudden" approximation which corresponds to the 
assumption that the electron acquires the photon energy in- 
stantaneously and emerges from the crystal surface immedi- 
ately. As a consequence, photoemission is associated with 
electrons near the crystal surface. In addition, only the mo- 
mentum component in parallel with the surface is conserved. 
It follows that ARPES is ideal for layered materials. The en- 
ergy levels involved in the ARPES process are shown in Fig- 
ure[Tt. Here, and throughout the paper, we define the quantity 
i?k coiTesponding to the dispersion of the paired fermions in 
terms of the usual BCS expression 



i?k= v/(ek-/x)2 + A2(T) 



(1) 



Because of the large photon energy hv, compared to the 
electron energy scale inside the crystal, the final state of the 
photo-emitted electron is essentially free so that the energy 
conservation constraint is given by Ei = Ej — hi', where 
Ef — k2/2TOe is measured with an energy analyzer. Here 
me denotes electron mass. In turn, the momentum (in the 
known direction) has magnitude k — yJ2raQEj. The ARPES 
spectrum is given by 



/P'^°*°(k,c.) = Afo(k,j.)A(k,c.)/H 



(2) 



where Afo(k, z^) is a matrix element which depends on the 
photon energy. Apart from the matrix element and the Fermi 
function /(w), one sees that ARPES measures the electronic 
spectral function. 

The energy levels involved in an RF transition are shown in 
FigHJ'- Here fi^ is the RF frequency for exciting a free atom 
from hyperfine level 2 (maroon line) to level 3 (green line). 
We neglect final state effects, which will be discussed later. 
A significant difference between an RF and ARPES transition 
is that in the RF case a dominantly large fraction (i7 1) of the 
photon energy is converted to excite a fermion from one in- 
ternal state to another. As a consequence, the excited atoms 
do not have a substantially higher kinetic energy so that they 
do not leave the bulk gas immediately after the transition un- 
til they are deliberately released. The energy zero for an RF 
transition is more conveniently chosen to be the bottom of the 
free atom band of state 2. In this convention, the final state 
energy is Ej = il^ + ek, where ek = k'^ /2m, and the initial 
state energy is Ei = ±i?k + A* for the two branches shown in 
Fig. lb. Therefore, the same energy conservation constraint 
emerges hv — Ej ~ Ei. Finally the RF current (which will 
be derived in Section HTCl l is 



/o^-^(k,H= MA(k,..)/H 

zvr 



(3) 



uj—e-\^ — 5v 



where |Ta; p is a tunneling matrix element and, the momentum 
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Figure 1: Figure on left: Energy levels in an ARPES transition. In a paired system there are two fermionic states which contribute to the 
photoemitted current. These correspond to the red and black curves. The upper branch (red curve) will not be occupied until the temperature is 
high. Here a tight binding dispersion tk is assume for the underlying non-paired initial state (magenta curve). The cyan dashed line indicates 
the Fermi level of the electrons and the green solid line the dispersion of outgoing electrons. Figure on right: Energy levels in an RF transition. 
f^L is the RF frequency for exciting a free atom from hyperfine level 2 (maroon line) to level 3 (green line). Q'j^ is the same energy but 
measured relative to the respective chemical potentials. The black and red curves are the dispersion of the particle and hole branch of a paired 
atom in level 2, with energy level given by =Fi5k + /i, respectively. 



integrated current, which is the more widely studied form, is 



^ 27r 

k 



A(k,c.)/(a;) 



(4) 



We note an important contrast with Eq. (O because here there 
is the restriction uj = e^-- Sv which (apart from the matrix 
element effects) serves to differentiate the photoemission and 
RF responses. 



B. Overview of the Literature on RF Experiments 

Experiments and theory have worked well hand in hand in 
developing an understanding of the so-called "RF pairing gap 
spectroscopy" in the atomic Fermi gases. This class of exper- 
iments was originally suggested by Torma and ZoUer |30] as 
a method for establishing the presence of superfluidity. In this 
context an equation equivalent to Eq ^ was derived. Later 
work Oil I32I1 . made the observation that these RF experi- 
ments, which reflect the spectral function A(k, lo), would ob- 
serve a pairing gap A(r) which may be unrelated to supercon- 
ducting order (except in the strict BCS regime). This was the 
beginning of a recognition that a pseudogap would be present, 
which is associated with stronger-than-BCS attractive inter- 
actions. Moreover, this pseudogap appears in the "fermionic 
regime", that is, when the fermionic chemical potential is pos- 
itive iiid . 

An experimental ground breaking paper |33] reported the 
first experimental implementation of this pairing gap spec- 
troscopy in ^Li over a range of fields corresponding to the 
BCS, BEC and unitary regimes. Accompanying this paper 



was a theoretical study 113411 by Torma and co-workers based 
on the BCS-BEC crossover approach introduced earlier ll3ll] . 
but, importantly, generalized to include trap effects. This the- 
oretical scheme is the one that will be the focus of the present 
paper The calculations showed reasonable agreement with 
experiment, and subsequent work lisll presented more quan- 
titative comparisons of the spectra along with theoretically- 
inferred estimates of the temperature, based on an adiabatic 
sweep thermometry ll36ll . Some of the first evidence that one 
was, indeed, observing a pairing gap (or pseudogap) in the 
normal phase was presented in Reference llJf], based on this 
same thermometric approach and the data of the Innsbruck 
group [33]. 

In an important contribution Yu and Baym pointed out fy% 
that the theoretical framework described above and summa- 
rized in Eq© missed what have now become known as "final 
state effects". Moreover, this could be seen most clearly in 
sum rule constraints on the RF spectra. These final state ef- 
fects can be understood as follows. Assume as the right panel 
of Figure [T] that the condensed phase involves pairing among 
hyperfine channels 1 and 2 and that the excited atomic state is 
associated with hyperfine level 3. While the attractive interac- 
tion gi2 drives the pairing, the excited atoms in 3 will also ex- 
perience a residual interaction 513, which may modify the RF 
spectra. In this way, these final state effects yield corrections 
to the lowest order current, shown in Eq (01). Interestingly, the 
sum rule, now known as the "clock shift" sum rule 1 38] shows 
that the first moment of the current sums to an internally con- 
sistent value, rather than a pre-determined constant. This will 
be discussed in Section lTlI Al 

A new set of groundbreaking experiments from MIT have 
introduced a powerful way of exploiting and enhancing RF 
spectroscopy first via tomographic techniques |39]. With the 
tomographic scans, the complication of studying the spectra 
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in a trapped configuration can now be removed, so that the 
system is effectively homogeneous. Also important was the 
demonstration that the entire collection of ^Li superfluids with 
hyperfine levels 1 and 2 paired, as well as 1 and 3 as well 
as 2 and 3, are stable and can be probed in RF spectroscopy 
with variable RF transitions, H. l (defined in the right panel of 
Figure [T]). In this way one has, in conjunction with a larger 
complex of superfluids, a way of tuning final state effects. 
Moreover, it was hoped that a proper choice of the superfluid 
and the RF transition can reduce the importance of these final 
state corrections and allow one to consider the simpler theory 
ofEq®. 

The theoretical challenge of incorporating final state con- 
tributions has become very topical, in large part because of 
the existence of data in effectively "homogeneous" systems 
through these tomographic techniques. It is only in the ab- 
sence of a trap that one can readily handle the higher order 
terms introduced by Yu and Baym iIstIi . With these correc- 
tions to Eq^ one may have a better opportunity to quanti- 
tatively fit the RF spectra. Very nice calculations ll40l l4lll of 
in the homogeneous case consider the T w superfluid 
and good agreement with experiment has been demonstrated 
ll40ll . Subsequent work ll42ll has addressed the entire range of 
temperatures where one can probe the RF contributions as- 
sociated with pre-existing thermally excited quasi-particles. 
These are shown as a second branch of RF transitions in the 
right panel of Figure [T] The body of work [35, 43] at general 
temperatures T makes the important point that the presence 
or absence of superfluid order (as long as T < T*) will not 
lead to fundamentally different physics. This observation is in 
contrast to alternative calculations ll40ll4Tll44ll which consider 
only the T « superfluid and/or separately the normal phase. 

Along with these new developments has been an experi- 
mental and theoretical focus on population imbalanced gases 
llH m m 113 m. The observation (4^ that extreme im- 
balance may drive the system to an exotic normal phase has 
captured the attention of the community. This exotic phase 
appears to be associated llsol IsTll with the binding of a small 
number of reverse spins to the majority states and this signa- 
ture is consistent with RF experiments, as shown theoretically 
ll44l l52ll. It should be stressed that this binding is not the same 
as pairing which is a macroscopic many body phenomenon. 
But it may, nevertheless, smoothly evolve into pairing as one 
varies the concentration of reverse spins L53J, and in this way 
diminishes the population imbalance. 

With the growing appreciation for final state effects, 
an interesting controversy has recently emerged concerning 
slightly different data obtained on the 12 superfluid at unitar- 
ity. This involves the original Innsbruck experiment and 
more recent data from the MIT group |24]. The latter series 
of studies have led the authors to inquire as to whether the 
pairing gap observations reported in Reference fy^ might in- 
stead be associated with final state effects. We comment on 
this possibility in Section |V CI of the paper, where we argue 
on behalf of the original interpretation in Reference ll33ll . 

Finally, recent experiments on "'"K from the JILA group 
ll20ll have now demonstrated that it is possible to measure the 
spectral functions directly using momentum resolved RF pair- 
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Figure 2: (Color online) Temperature dependence of the excita- 
tion gap from the ARPES measurement for optimally doped (filled 
black circles, Tc = S7K), underdoped (red squares, Tc — S3K) 
and highly underdoped (blue inverted triangles, Tc = lOA') single- 
crystal BSCCO samples (taken from Ref. 1 54]). There exists a pseu- 
dogap phase above Tc in the underdoped regime. 

ing gap spectroscopy over a range of magnetic fields through- 
out the BCS-BEC crossover. In these recent experiments 12(1 
the momentum of state 3 atoms is obtained using time-of- 
flight imaging, in conjunction with 3D distribution reconstruc- 
tion. Since the 3D gas in a single trap is isotropic, detailed 
angular information is irrelevant. There is a substantial ad- 
vantage of using K'^'^ over the more widely studied ^Li since, 
for the usual Feshbach resonance around 202 G, there are no 
nearby competing resonances to introduce complications from 
final state interactions ][3l,|4^|4ll|4l|5l. This powerful tool, 
which we have seen has a strong analogy with ARPES spec- 
tra, opens the door for testing the fundamentals of the many 
body theory which underly this BCS-BEC crossover. As we 
show later in Section IV CI it also helps to remove ambiguity 
plaguing the interpretation ll24l [33|] of momentum integrated 
RF experiments by establishing a clear dispersion signature of 
pairing. 

C. Key Features of ARPES Data on Cuprates 

We outlined earlier three issues around which much of the 
cuprate photoemission studies can be organized. These are 
characterizations and modelling of the fermionic self energy, 
of the pseudogap phase and of the effects of coherence as the 
superconductor passes from above to below Tc- Figure |2] is 
a plot showing the behavior of the excitation gap which ad- 
dresses the first of these three issues. Plotted here is the pair- 
ing gap inferred from the leading edge in the photoemission 
experiments as a function of temperature. The temperature T* 
can be read off as the temperature where the gap first appears. 
The three different curves correspond to three different doping 
concentrations which one can interpret in the framework of 
BCS-BEC crossover as corresponding to three different val- 
ues for the attraction strength, since they correspond to three 
different values for the pairing onset T*. 
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Figure 3: (Color online) Temperature dependent photoemission 
spectra from optimally doped Bi2212 (Tc — 91K), angle integrated 
over a narrow cut at (7r,0). Inset: superconducting peak intensity 
vs temperature. After Ref.|21]. Note the sharpening of the peaks as 
temperature is lowered below Tc- 



Several key points can be made. The transition tempera- 
tures for phase coherent order are not evident when one stud- 
ies the pairing gap, as shown in the figure except, perhaps, in 
the sample with the lowest T* corresponding to the highest 
doping concentration. In this case T* « Tc as in the BCS 
limit. In general the higher is T* the lower is Tc which is, in 
fact, consistent with what one expects for BCS-BEC crossover 
on a lattice llssi |56[ IstIi . Because the behavior varies so 
smoothly from above to below Tc one says that the pseudo- 
gap smoothly evolves into the excitation gap of the superfluid 
phase (sometimes called the "superconducting gap"). Indeed, 
independent experiments show that these two have the same 
(i-wave symmetry. 

Ideally, one would like to obtain the analogous plots which 
show the temperature dependence of the pairing gap in the 
cold Fermi gases, using RF spectroscopy in the BCS, unitary 
and BEC regimes. It should be noted that none of the cuprate 
curves represent the BEC case. Even thou gh t he ratio T* /Tc 
can be quite large due to lattice effects ll55ll58ll . the high tem- 
perature superconductors are all in the fermionic regime. 

In Figure[3]the cuprate photoemission spectra near optimal 
doping (Tc ~ 91 K) are plotted for a range of different tem- 
peratures in order to exhibit the effects of emerging phase co- 
herence. This figure represents angle integrated spectra over 
a cut near the d-wave anti-node (where the gap is largest). 
What is most striking here is the fact that a sharp quasi-particle 
peak emerges only below Tc- Above Tc there is a pairing gap, 
but it is associated with a relatively poorly defined (gapped) 
quasi-particle. In the BCS-BEC crossover scenario, one might 
view this as representing short lived, but non-condensed pairs, 
which only below Tc can become long lived and stable. 

Also of interest in the figure is a feature known as the peak- 



dip-hump structure which is associated with superconducting 
coherence. There is still some controversy over the origin of 
these effects, but some ||23t] have correlated them with specific 
bosonic modes which couple to the fermions and appear in the 
self energy. 

Again, it would be interesting to have more complete anal- 
ogous studies using RF spectra on the cold gases as the sys- 
tem varies from above to below Tc- Just what are the pre- 
cise signatures of superfluid coherence and is there evidence 
that short lived non-condensed pair states become longer lived 
below Tc, needs to be addressed. Here one must (perhaps 
through tomography) overcome the complexity introduced be- 
cause these gases are contained in a trap. 

II. GENERAL THEORETICAL BACKGROUND 
A. BCS Leggett T-matrix Theory 

This paper will address the theory behind RF spectra and 
photoemission in the cuprates in the context of one particular 
approach to BCS-BEC crossover based on the BCS-Leggett 
ground state. Here, however, we generalize to finite temper- 
atures T. There is an alternative approach (aO, '4?! based on 
the Nozieres Schmitt-Rink scheme |58] and which involves 
another ground state. Because of the flexibility of the BCS- 
Leggett scheme which can readily be generalized to include 
trap effects within the superfluid phase, as well as popula- 
tion imbalance, we choose this alternative. Another major 
advantage (from our perspective) is that it is not plagued by 
issues associated with a first order transition lf59ll at Tc- These 
effects are related to analogous behavior in mean field theo- 
ries of the Bose gas. Interpretation of the cuprate data, which 
shows a smooth evolution through Tc, would be problematic 
in the presence of first order effects. 

We briefly summarize the key equations which emerge 
from our T matrix scheme. Within the present approach there 
there are two contributions to the full T-matrix t = tpg + tgc 

where tsc{Q) = — j^S{Q), where Age is the superfluid (sc) 
order parameter Similarly, we have two terms for the fermion 
self energy ^{K) = + ^pgiK) = Y.q t{Q)Go{Q - 

K). Here K and Q are four- vectors. It follows then that 

Esc(k, cj) = cj+eil-jj • Throughout this paper the label p5 
corresponds to the "pseudogap" and the corresponding non 
condensed pair propagator is given by 

tpg[Q) = Ul[l + Ux{Q)l (5) 

where the pair susceptibility x{Q) has to be properly chosen 
to arrive at the BCS-Leggett ground state equations. We im- 
pose the natural condition that below Tc there is a vanishing 
chemical potential for the non-condensed pairs 

^^pair = (6) 

which means that tpg{Q) diverges at Q = when T < Tc- 
Thus, we approximate [60, 61] Y.pg{K) to yield 

i:pg{K)^-Go{-K)Al^ T<Tc, (7) 



6 



with 



B. Generalization To Include Population Imbalance 



(8) 



It follows that we have the usual BCS-like form for the self 
energy E(k, uj) « A^/[uj + ek - m] with T <Tc with 



A2(T) = AUT) + AI{T). 



(9) 



As is consistent with the standard ground state constraints, 
Apg vanishes at T = 0, where all pairs are condensed. 

Using this self energy, one determines G and thereby can 
evaluate tpg. Then the condition that the non-condensed pairs 
have a gapless excitation spectrum (fipair = 0) becomes the 
usual BCS gap equation, except that it is the excitation gap A 
and not the order parameter Age which appears here. We then 
have from Eq. (|6]l 



For consistency we must take for the pair susceptibility 



x{Q)^yGoiQ~K)G{K). 



K 



(10) 



(11) 



Here G = {Gq^ — E)"^ and Gq are the full and bare Green's 
functions respectively. 
Similarly, using 



K 



one derives 



E 



, £k - M 
Sk 



(12) 



(13) 



which is the natural generalization of the BCS number equa- 
tion. The final set of equations which must be solved is rather 
simple and given by Eq. ([8]), Eq. ( fTOl ), and Eq. ( fT3] ). Note 
that in the normal state (where iJ^pair is nonzero), Eq. O is 
no longer a good approximation, although a natural extension 
can be readily written down |43]. 

We stress that the approximation in Eq. O is not central to 
the physics, but it does greatly simplify the numerical analy- 
sis. One can see that correlations which do not involve pairing, 
such as Hartree terms are not included here. This is what is 
required to arrive at the BCS-Leggett ground state. It should 
be clear that, in principle, the T-matrix approach discussed 
here is more general and that in order to address experiments 
at a more quantitative level it will be necessary to go be- 
yond Eq. O. Various groups |40, 44] have included these 
so-called "CqGo" contributions to the pair susceptility xiQ) 
in the T-matrix. These are particularly important when the 
pairing is weak at high temperatures, or for imbalanced gases 
II44II51II52II . and, like the Hartree corrections, will ultimately 
have to be included. 



One major advantage of the BCS-Leggett approach is that 
it is straightforwardly generalized to include a population im- 
balanced superfluid. We begin by summarizing the general 
equations associated with the so-called "Sarma state" which 
corresponds to a uniformly polarized BCS superfluid. This is 
to be distinguished from the phase separated state ll45ll . 

The gap equation is given by 







1 V- 

77 + E 



U ^ 2Ek 

k 

Here we define the average 

f(x) = [f{x + h) + f{x-h)]/2, 



(14) 



(15) 



where f{x) is the Fermi distribution function. In addition 
we define ^ = (/Lt| + and h = {^-^ — fJ.i)/2, E\s_ = 

Ve^TA^, Ey,^ = -h 
^k = Ck - M- 



i?k and i?ki = h + _Ek, where 
There are now two number equations given by 



Ey 



-fiEu 



6n = J2ifiEk-h)-f{E^ + h)], 



(16a) 



(16b) 



where n = n| +n| is the total atomic density, Sn = n^ — ni > 
is the number difference and S = 5n/n is the polarization. 
Here the coefficients u^, — (1 ± ^k/^'k)/2 are formally 
the same for both the polarized and unpolarized systems. 



Finally, one has to recompute A: 



pg 



Al^{T) = A\T) - Al{T) = -J2 tM 



(17) 



This quantity can be obtained directly from an equation in 
the same form as Eq. ^ except that the pair susceptibility 
appearing here satisfies 



xiQ) = ^[xniQ) + xi^iQ)] 



(18) 



where as before we have the product of one dressed and one 
bare Green's function 



xniQ) = Y.^oi{Q^K)GUK) 

K 

XiT(Q) = Y.^oi{Q - K)G^{K) 

K 



(19a) 
(19b) 



Further details are presented in Reference II43I1 . 



C. Linear Response Theory and RF 

In the RF experiments [62], one focuses on three different 
atomic hyperfine states of the ^Li atom. The two lowest states. 
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|1) and 1 2), participate in the superfluid pairing. These corre- 
spond to t and j. The higher state, |3), is effectively a free 
atom excitation level; it is unoccupied initially. An RF laser 
field, at sufficiently large frequency, will drive atoms from 
state |2) to |3). 

We begin with the usual grand canonical Hamiltonian H — 
liN which describes states |1) and |2). We have described in 
Section HTAl the procedure for handling pairing correlations in 
this 12 channel. The Hamiltonian describing state |3) is given 
by 



the spectral function associated with the superfluid component 
1 2). Then the lowest order RF current can be written as 

= E ^(k, - i') [/(Ck -1^)- /(ek,3)] , 

k 

(21) 

In practice, state 3 is unoccupied, thus the second Fermi 
function in brackets vanishes. In this way we have for the 
momentum-resolved current 



/o(k,i.) = -— A(k,ek-j^)/(ek-i^). (22) 

ZTT 



where eu is the atomic kinetic energy, Cg is the annihilation 
operator for state |3), and /ig is the chemical potential of |3). 

In addition, there is a transfer matrix element Tk.p from |2) 
to 1 3) given by 



'^3,p'^2,k 



h.C.) 



k,p 



For plane wave states, Tk.p = 723'^(qL + k — p)(5(wkp — 
fii). Here ql ^ and JIl are the momentum and energy of 
the RF laser field, and Wkp is the energy difference between 
the initial and final states. In what follows we will set the 
magnitude of the tunneling matrix to unity, without loss of 
generality. It should be stressed that unlike conventional SN 
tunneling, here one requires not only conservation of energy 
but also conservation of momentum. 
The RF current is defined as 

I = -{N2)=-t{[H,N2]). 
Using standard linear response theory one finds 

/H = --Im[i?^(17L+M-M3)]. 

TT 

Here we introduce the retarded response function D^{uj) = 
D{iuJn LO + iO+). 

At the lowest order of approximation the linear response 
kernel D can be expressed in terms of single particle Green's 
functions as 

K 

where K = (k, a;„) and Q = (0, ri„). (We use the convention 
h = ks = I)- The Green's function can then be expressed in 
terms of spectral functions. After Matsubara summation we 
obtain 



^ /'de^A(k,e)A3(k, £-)[/(£-) 



fie)] , (20) 



with e = e + i^ + /i — fi^, v is the RF detuning and J{x) is the 
Fermi distribution function. 

The spectral function for state |3) is v43(k, e) — 27r(5(e — 
^k.3)- We see from the above equation that these RF exper- 
iments depend on A(k, e) = —2 ImG(k, e + iO+) which is 



D. BCS-Leggett Model for Self Energy 

The current /o(i') at the leading order level depends on the 
fermionic spectral function, which, in turn depends on the 
fermionic self energy. In this section we discuss the nature 
of the self energy which will enter into an analysis of both 
photoemission and RF spectroscopy. 

To arrive at the BCS-Leggett ground state equations, we 
have seen that the self energy is given by I](k, lj) = 
S.,^(k, w) 



S(k,c^) 



S]pg(k, w) where 

k,sc 



w + ek — M 

k.sc 



Spg(k,w) 



(23) 



— ,T<T^ (24) 



A£ 



a; + ek - /i 



(25) 



These equations follow, provided one makes the approxi- 
mation contained in Eq. (I?). In invoking this approximation 
we are in effect ignoring the difference between condensed 
and non-condensed pairs which cannot be strictly correct. The 
simplest correction to Spg (which should apply above and be- 
low Tc) is to write 



Spg(k,w) 



A^ 



+ Ek - + «7 



So(k,u;) 



(26) 



Here the broadening 7 7^ and "incoherent" background con- 
tribution Eo reflect the fact that noncondensed pairs do not 
lead to true off-diagonal long-range order While we can think 
of 7 as a phenomenological parameter in the spirit of the high 
Tc literature ll63l l64ll there is a microscopic basis for consid- 
ering this broadened BCS form ll60ll65ll . The precise value of 
7, and even its T-dependence is not particularly important for 
the present purposes, as long as it is non-zero at finite T. For 
simplicity we will take 7 as a temperature independent con- 
stant throughout this paper By contrast Ssc is associated with 
long-lived condensed Cooper pairs, and is similar to Spg but 
without the broadening. 

It is important to stress that this same self energy model has 
been applied to describe the spectral function in the pseudogap 
116311641 l66tl and the superfluid phases |67] of the high temper- 
ature superconductors, where here So(k, oj) = —iVq is taken 
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to be an imaginary constant. In the cuprate literature, it has 
been argued that Eq. ( |26l ) is appropriate to the normal phase 
and the onset of coherence coincides with a dramatic decrease 
in 7 below Tc- This leads to a subtle and important contro- 
versy about which we wish to comment more in the context 
of Section [Vll Our own perspective is that Eq. (|23] | in con- 
junction with Eq. ( |26] l is the appropriate starting point. That 
is, there are two contributions to the self energy below Tc and 
only one above. Thus one should not argue that 7 precisely 
vanishes at T^. but rather there is a continuous conversion from 
non-condensed to condensed pairs as T is lowered within the 
superfluid phase. The non-condensed pairs below Tc have fi- 
nite lifetime while the condensed pairs do not. 

The resulting spectral function, based on Eq. ( |26] l and 
Eq. ( I23] ) is given by 



^(k,e) 



2A2 7(e + ^k)2 



(e + a)^(e2-£;2)2+y(e2_^2_^2j2 



(27) 

Here, for convenience we do not show the effects of the So 
term. Above Tc, Eq. (|27| i can be used with A^c = 0. It can 
be seen that this spectral function at all k contains a zero at 
e = — below Tc, whereas it has no zero above Tc- This 
means that a clear signature of phase coherence is present, as 
long as 7 7^ 0. 

This analysis may be readily generalized to include the ef- 
fects of population imbalance. We have for the spectral func- 
tion of the minority 



^2(k,e) = — 



2A^^7(^'+a)^ 



with e' = e — h. 

Similarly, the spectral function for the majority is 



(28) 



Ai(k,e) 



2Ag,7(," + a) 



(e" + Ck)2(e"2 - Elf + 72(e"2 - - A^y 



(29) 

with e" = e + h. 

There are instances where it is problematic to include the 
effects associated with the finite lifetime 7. This occurs when 
we compute the effects of final state interactions. At this strict 
"mean field " (mf) level we drop the factor 7, thereby, losing 
the distinction between condensed and non-condensed pairs. 
In this case the spectral function (which we display here for 
the polarized gas case) is associated with majority (1) and mi- 
nority (2) contributions: 



A^f{k, iv) = 2n[ulS{Lu - Sk.i) + vlSiuj + iJk.t)], 

with i?k,| = Ei^ — h and E'k.x — Ey^ + h. We stress 
again, however, that here A contains non-condensed pair ef- 
fects through Eq. (|9|l. 



E. Analytical Results for the Leading Order RF Current 

It is possible to obtain analytical results for the leading or- 
der current at general temperatures T in this strict mean field 
theory. Here one can integrate over the momentum to find 



1 A_2 

^2 

^^f{Eo)ko, (z.2<i^<0) (31) 



^oM = ^ — [^-fiEo)]ko, i'y>i^i) (30) 



with 



Eo 



2v 



^5 = M + 



2v 



The frequency regimes associated with the negative and pos- 
itive continua are given by -~{\/ n"^ + A^ + ^j) < v < Q 
and V > ^ [P- + A^ — /i. In the above equations V2 = 

-(Vm^ + A2 + ^) and vx = \J ii^ + A^ - ^. It can be 
seen that there are contributions for both negative and posi- 
tive detuning. At strictly zero temperature, the Fermi function 
vanishes, and we have only the positive continuum 



1 A^ / _ A2 



(32) 



It is possible to write down a closed form expression for the 
polarized case at general T as well. We do not show it here 
because the expressions are more cumbersome. We defer this 
to AppendixlAl 



F. Behavior in Traps In Absence of Final State Effects 

Once the trap is incorporated, one has to solve for the cur- 
rent at each position r and then integrate in the form 



(33) 



where here na-{r) represents the particle density within the 
trap and a = 1 , 2 are the different hyperfine levels of the 
superfluid. To handle trap effects we assume a spherically 
symmetrical harmonic oscillator potential V{r) — mOj^r'^/2. 
The density, excitation gap and chemical potential which vary 
along the radius can be determined I143il using the local density 
approximation (LDA). 

It should be stressed that the density and gap profiles (ricr (r) 
and A(r)) in general involve pseudogap or non-condensed 
pair effects. The strict mean field theory, which often gives 
a reasonable approximation to the spectral functions, is not 
adequate for obtaining these trap profiles. Thus even when 
analyzing tomographic RF data one has to include the full ef- 



fects of these pair excitations Il43ll . effectively through Apg{r) 
and non-zero ^JLpair{'^')■ 
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K, 2 ^12 




Figure 4: Feynman diagrams for the RF response function D{Q). The left bubble is the lowest order Do, whereas the right diagram, Dal, is 
associated with final state effects. Here thin (thick) lines stand for bare (full) fermion propagators, the dashed line for ti2, approximated as the 
condensate, and double wiggly line for tri- The numbers in blue indicate the hyperfine levels. Here Q = (in„, 0) for the RF field. 



III. FINAL STATE EFFECTS IN HOMOGENEOUS 
UNPOLARIZED SYSTEM 

We now turn to the inclusion of final state effects which go 
beyond the leading order diagram. It is complicated to han- 
dle these contributions for the inhomogeneous case. Thus we 
focus here on treating the RF current as if the system were 
homogeneous. We formulate the finite T, RF problem using 
a diagrammatic scheme. The full diagram set for the RF re- 
sponse function, D{Q), is shown in Fig. 1. The leading order 
term Dq appears as the first term on the right hand size and the 
second contribution is associated with the Aslamazov-Larkin 
(AL) diagram (called Dal)- The full RF current, given by 
the retarded response function, is = — (l/7r)ImZ)^(il), 
where il = i/ + fi — 

The approximation compatible with Eq. (|7]) is effectively 
equivalent to treating the Dal in Fig. 1 at the BCS mean-field 
level, leading to the opposite momenta ±K for particles 1 and 
2 in the diagram. Dal{Q), depends on A, not A^c, and in- 
corporates final-state effects via the interactions gi2 between 

1 and 2 and 513 between 1 and 3. We neglect the effects aris- 
ing from the interaction between 2 and 3. This is consistent 
with the approach in Ref. [38]. This second term has appeared 
previously in studies of the superfluid density ll68ll . Our for- 
mulation of the finite T, RF problem can be made compatible 
with the diagrams in Ref. 1,40.1 . although attention in that paper 
was restricted to very low temperatures. Our diagrammatic 
scheme reduces at T = to the approach of Ref. |41]. 

In order to evaluate the AL term, we begin by writing out 
the relevant T-matrices 

tiiiQ) = gr2' + E^i(^)^2(Q-^) (34) 

K 

tu\Q) = g^i + Y.^i{K)Gl{Q ~ K) (35) 

K 

where 512 and 313 parameterize the interaction between 1 and 

2 and 1 and 3, respectively. We can also introduce the s-wave 
scattering lengths, ais (and 012) in the 1-3 (and 1-2) channels, 
respectively. 

Thus 



where 



Ana 
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xi3(Q) 



(36) 



1 ~ fjEu) - /(ek,3) 

1 

ifln — E]^ — ^k,3 



/(i?k) - /(Ck,3) 



iQ.n + Els, — ^k,3 

The AL diagram yields 

Dal{Q) = [Y,FiK)Gt{K + Q) 

K 

where 

F{K) = -AG'^^\K)Gi^\-K) = 



1 



(37) 



tMQ). (38) 



A 



This contribution can be rewritten 

Dal{Q) = DliQ)h3{Q). 

where we have defined 

D2{Q)^j:KFiK)GUl< + Q) 



(39) 



(40) 



E 

K 



A rl-/(£;fe)-/(6,3) f{Ek)- f{ik,3) 



2Ek 



iQ.n - Ek - i 



k,3 



iVLn + Ek - i 



k,3 



(41) 



Then the full set of diagrams shown in Figure|4]can be com- 
bined to yield 



DiQ) = Dq{Q) 



(42) 



After analytical continuation and change of variables, we 
have 57 ± E]^ — ^j. 3 = v ± E\^ — Importantly, the de- 
nominators here are the same as those which appear in ti2. 
Furthermore, at = 0, /(Cfc.s) is cancelled out so that 

^ (0) = (51-3' - 9i2 )+ti2i0)^ 9ii - 9i2 ■ (43) 

It follows that the complex functions Do{Q), Xi.3(Q), 
D2{Q) are the same as their wave function calculation coun- 
terparts |41] when the pairing gap A is chosen to be order pa- 
rameter Age and T = 0. It is not fl that should be identified 
with the experimental RF detuning. 
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After some straightforward algebra (with details in Ap- 
pendix A), we find for the RF current 





1 


1 ' 






.512 


513. 




1 


f 


m 


m \ 






V47rai3 


47rai2 / 






1 


m 


m 


2^2 

— Imtf3(;. 




TT 


47rai3 


47rai2 



Moreover, in the special case, when 013 = ai2 then — 
(n2 — n3)6{h'). Equations (l44l i are the central result. It should 
be clear that final state effects in the RF current directly reflect 
the T-matrix in the 1-3 channel. In general, features in the RF 
spectra derive from the poles and imaginary parts of Dq{Q), 

XiAQl andi^2(g). 

The spectrum may contain a bound state associated with 
poles at vq in ^13, as determined by t^-l{vQ) ~ 0. This leads 
to the so called "bound-bound" transition. In addition, there is 
a continuum associated with both the numerator and denomi- 
nator in the first of Eqs. (|44] |. with each contribution spanned 
by the limits of = ^k±Ek, i.e., -{sj + A2 + ^) < 1/ < 
and V > y^fj^ + A^ — /i. The continuum at positive frequen- 
cies is primarily associated with breaking a pair and promot- 
ing the state 2 to state 3. This represents the so-called "bound- 
free" transition. On the negative detuning side, the continuum 
is primarily associated with promoting to state 3 an already 
existing thermally excited 2 particle. The spectral weight of 
the negative continuum vanishes exponentially at low T. 



A. Sum Rules 

Of importance, in assessing any theoretical framework for 
computing the RF current are the two sum rules associated 
with the total integrated current and the first moment or "clock 
shift" [38]. Using the Kramers-Kjonig relations between 
Retfg and Imi{|, we prove in Appendix ID] that, not only in 
the ground state, but also at finite temperature, Eq. (l44l l satis- 
fies 



dv I{v) = 712 - , 
dvvl{v) = 



m 
47r 



1 

ai2 



1 

ai3 



(45) 
(46) 



where 712 and n-i{— 0) are the density of state 2 and 3 atoms, 
respectively. In this way we find for the clock shift a result 
which we write (for general polarizations, associated with the 
subscript a) in the form: 



/ dvlc{v) - ns At: 



1 

ai2 



1 

ai3 



(47) 



In the unpolarized case, this agrees with Ref. 113811 . This 
sum rule is satisfied only when ai3 7^ and when both di- 
agrammatic contributions are included. It is easy to show 




Pseudogap 



Superfluid 



■ stronger attraction 



Figure 5: Phase diagram showing Tc and T* for homogeneous s- 
wave Fermi gas superfluid (left) and for d wave superfluid on a quasi 
two dimensional lattice. From Reference |69]. Note that the BEC 
asymptote is finite in a Fermi gas and zero in the lattice case. Because 
the lattice phase diagram shows similarity to that of the cuprates, in 
future experiments it will be important to study the d-wave general- 
ization of the attractive Hubbard model on an optical lattice. 
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Figure 6: This summarizes the phase diagrams for polarized gases 
with and without a trap and with and without pairing fluctuations. 
The figure is based on Refs. jksll and (4^. The figures on the left are 
for homogeneous case and on the right for trapped case. The mean 
field figures at the top show the reported tri-critical point. However, 
fluctuations (in the lower two plots) depress the superfluid phases. 



that at large IqIv) ~ v^'^/'^, Imtp_^~ so that 

^ v^^/"^, in agreement with Ref. ll40ll . Clearly, the first 
moment of I{h') is integrable, whereas the first moment of 
Iq{v) is not. Finally, Eq. (|44] | reveals that the spectral weight 
(including possible bound states) away from v — Q will dis- 
appear when the gap A vanishes. 



rV. PHYSICAL PICTURE AND IMPLICATIONS 

In this section we lay the groundwork for a comparison be- 
tween theory and experiment, which is presented in the fol- 
lowing section. We address the various phase diagrams for 
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the population balanced Fermi gases, including the (d-wave) 
lattice case, as well as for the imbalanced systems. We an- 
alyze a pedagogically useful set of figures which lay out the 
general behavior of the RF spectra with and without final state 
effects and with and without a trap. Importantly we compare 
photoemission-based plots for the same parameter set as RF- 
based plots and address the key signatures of emerging super- 
fluid coherence as one goes from above to below Tc. 



A. Phase diagrams 

The relevant phase diagrams to be used and referred to in 
our RF calculations have been obtained elsewhere. Shown 
here are the curves for Tc and T*. The latter represents the 
boundary between the pseudogap and normal phases. Fig- 
ure|5]compares the phase diagram for an s-wave paired Fermi 
gas (left) and for a d-wave paired fermion system (right) on a 
quasi two dimensional lattice. The s-wave gas case is closely 
analogous to results obtained using the approach of Refs. |d8 ] 
and ILZOD. The d-wave case was discussed earlier in Ref. [TlJ] 
and more recently in the context of optical lattice calculations 
in Refs. ifssll and ll69ll . The seminal Nozieres Schmitt-Rink 
paper pointed out a key fact which identifies a notable differ- 
ence between the lattice and gas cases: the BEC limit has an 
asymptote of Tc ^ in the case of a lattice, whereas it is finite 
in a gas. Thus there is a relatively larger separation between 
T* and Tc when fermions are present on a lattice as shown 
in the figure. We note that the d-wave case has a number of 
features in common with the counterpart phase diagram of the 
cuprates |55, 69, 71]. 

The phase diagrams for polarized (s-wave) Fermi gases are 
shown in Figure |6] The four panels correspond to the effects 
of including (or not) a trap and to the effects of including 
(or not) pairing fluctuations beyond strict mean field theory, 
which enter in the theory through the parameter A^^. When 
we discuss the RF behavior of polarized gases we will use the 
full beyond-mean-field theory phase diagrams, although some 
of the calculations of the spectral function are performed at the 
strict mean field level. 

Beyond the normal phase, there are three phases which ap- 
pear |43, 48]: the Sarma phase, a phase separated (PS) state 
and a pseudogapped normal state, as indicated. We note that 
the treatment of the normal component of the phase separated 
state does not include correlations beyond those accounted for 
by A. As a result, these calculations overestimate the range of 
stability of phase separation. This issue has been nicely dis- 
cussed in the theoretical literature ifsol IsTIl with implications 
for RF spectra as well [|44i 15211 . The Sarma phase should be 
considered as the more correctly treated here and one notes 
an important finding: that in the absence of a trap the regime 
of stability of the Sarma state is greated reduced. This more 
restricted stability (seen by comparing the two lower figures) 
is associated with the fact that the excess majority fermions 
can be accomodated more readily in different spatial regions 
in a trap. The maximum polarization of this homogeneous 
Sarma phase is around S — 0.2 which is close to that reported 
experimentally [i72i] . 
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Figure 7: This figure shows the behavior of the excitation gap as a 
function of temperature |73] for the Fermi gases at three different 
scattering lengths. This should be compared with Figure |2] for the 
cuprates. Arrows indicate locations of Tc. 
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Figure 8: This is a photoemission-like plot for a homogeneous uni- 
tary Fermi gas based on Eqs. ([26} and Here Tc = 0.27 and 
T* ~ O.STf. The figure shows that the onset of superfluid coher- 
ence leads to a sharpening of the peak structure. We take 7 to be 
0.5i?F. This figure can be compared with Figure|2] 



B. Comparison with the Cuprates 

We present in Figure [T] a plot of the excitation gap in the 
cold gases for three different values of the s-wave scattering 
length in units of 1/kpa which are near unitarity (a = 00) 
and on both the BCS and BEC sides. This figure should be 
compared with Figure [2] for the cuprates. Here the excitation 
gap is estimated using Eq. ( fTOl ) for all temperatures |73]. Also 
indicated on the curves is the value of the transition tempera- 
ture. This figure makes it clear that pseudogap effects, which 
are essentially absent on the BCS side of resonance, are very 
apparent at unitarity, where the Fermi gas has a positive chem- 
ical potential. In both the unitary and BEC cases, A is roughly 
temperature independent below Tc- 

Figure [8] represents a photoemission-Uke study, but for the 
parameters associated with a unitary (homogeneous) Fermi 
gas. Here the vertical axis plots the k integral of /p'loto^k, cu) 
based on Eq. (|2]l assuming a structureless matrix element Mq. 
This figure should be compared with the cuprate data in Fig- 
ure [3] The various curves correspond to different tempera- 
tures as indicated with Tc ^ 0.25Tf and T* w O.bTp. The 
self energy is based on Eq. ( |26] | for the non-condensed pair 
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Figure 9: Pedagogical figure siiowing typical RF spectra of unitary 
homogeneous gas at temperature somewhat below Tc- The various 
characteristic energy scales are labelled. Upper panel corresponds to 
absence of final state effects while lower panel includes final state 
effects in rather extreme limit of a (weakly) bound state in the neg- 
ative continuum. This is optimal situation for using the sum rules to 
extract A. 
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Figure 10: Based on previous figure, plotted here are pairing gaps 
A{T)/Ef vs temperature as computed exactly and as computed 
from a limited integration, using the sum rule. 
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Figure 11: This is the RF counterpart of Figure [8] for homogeneous 
unitary gas. The main body of the figure plots the higher T behavior 
and the inset shows the results at lower T when superfluid coherence 
is well established. Just as in photoemission, there appears to be a 
signature of this coherence in the RF spectra which is associated with 
a rather sharp threshold behavior, as seen in the inset. 
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Figure 12: Evidence that Tc may show up as subtle feature in RF 
spectroscopy. Based on the results of Figure [TT] plotted here is the 
RF current magnitude at zero detuning (for a homogeneous balanced 
unitary gas without final state effects) as a function of temperature. 



the superfluid phase (if one compares with the same value of 
T*). 



C. Overview of RF spectra: Homogeneous Case 



component with 7 — 0.25Ep and jEq = Fq = O.lEp- What 
is most notable about this figure is the progressive sharpen- 
ing of the "photoemission" peaks associated with the growth 
of coherence as T decreases. This same effect is seen in the 
cuprate data (Figure O. One notes here, however, that there 
is some shift of this peak position reflecting an increase in A 
with decreasing T, which is not seen in the cuprate data. This 
effect can be attributed to the fact the there is a substantially 
larger separation I55l [57ll69tl between T* and Tc in the case 
of a lattice (away from the BCS regime) than for a gas such 
as shown here. This is apparent in Figure |5] Thus, there is 
more temperature dependence found in the excitation gap of 



The top panel of Figure |9] represents a plot of character- 
istic RF spectra for a unitary gas without final state effects 
and at a moderate temperature below Tc- Here we use strict 
mean field theory. Indicated in the figure are the various en- 
ergy scales showing the location of the pairing gap A as well 
as the thresholds associated with the negative and positive 
continua. Spanned by ~{\J + A^ + /.t) < < and 
V > + A^ — fi these are indicated as V2 and vi. One 
can see that there is a substantial separation between the pair- 
ing gap value and the threshold . Thus, there is very little in 
the figure to suggest a way of extracting the pairing gap. This 
has presented a dilemma for the field. 
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One way to address this issue is to exploit the sum rule 
in Eq. ( |47] l which is appropriate provide one includes final 
state effects. In the lower panel we show the same spectra 
when final state effects are included. We have chosen a very 
special case for illustrative purposes in which a (meta-stable) 
bound state overlaps the negative continuum. This represents 
the most ideal example for exploiting sum rule constraints to 
extract the pairing gap. One can see here that because the 
bound state is in the negative continuum, the bulk of the spec- 
tral weight is confined to a narrow frequency weight spanning 
from to 0. 

In Figure [TO] we show the estimated values for the pairing 
gap A of a unitary gas obtained from the sum rule as inte- 
grated from V — ~2Ep io v = +2Ep compared with the 
exact pairing gap. The accuracy is within 10 %. To arrive at 
a case where the final state is on the BEC side of resonance 
is reasonably straightforward and the 13 superfluid, which ex- 
hibits this, is now well studied by the MIT group [24]. How- 
ever, we point out that for this 13 superfluid and for typical 
values of kp the bound state is deep and well removed from 
the continuum. By contrast, the case shown here results from 
a situation in which fc^? is increased from the currently quoted 
experimental values by about a factor of 10. While this may 
not be easy to achieve in the near future, it does point to the 
advantage of exploiting final state effects to focus the spectral 
weight in the more well confined, negative i/ regime. 

We plot in Figure [TT] the homogeneous spectra in the ab- 
sence of final state effects but now for the case in which we 
go beyond strict mean field theory and, thereby, differenti- 
ate the condensed and non condensed pairs on the basis of 
Eq. ( |26] |. We have chosen the same parameters as in Figure|8] 
In contrast to this earlier photoemission-based calculation one 
applies Eq. (|4]i to describe the RF experiments. Each of the 
curves corresponds to the same temperatures as their coun- 
terparts in Figure |8] It can be seen that the shape of the RF 
spectra is very different from that of photoemission (even in 
the absence of final state effects), because of the constraint in- 
troduced hy CO = — Siy which appears in Eq. (|4]l. There 
are two peak structures at higher T, even in this homogeneous 
situation, with the lower peak corresponding to the negative 
continuum. The temperature regime is separated into the su- 
perfluid phase in the inset and the normal phase in the main 
body. The differences between the two sets of curves supports 
the fact that there are signatures of superfluid coherence in RF 
spectroscopy. While, there are no abrupt changes at Tc, nev- 
ertheless, like the quasi-particle peak sharpening seen in the 
cuprates (in Figure |3]l one sees a distinctive sharpening of the 
positive I' threshold as the temperature crosses Tc- 

To quantify this. Figure [12] shows a plot of the zero detun- 
ing current for this balanced gas as a function of temperature. 
There is a feature at Tc which may make it possible, in princi- 
ple, to extract this transition temperature from high resolution 
RF experiments. This feature reflects the fact that there are 
subtle but distinctive features in the RF spectrum associated 
with superfluid coherence, as seen in Figure [12] Indeed, one 
infers from this previous figure that the negative continuum 
only begins to be perceptible as the temperature approaches 




0.5 1 1.5 2 



Figure 13: Contour plots of momentum resolved RF spectra for the 
homogeneous case at unitarity at T/Tc ~ 1.9. The top and bottom 
lines in the contour plots correspond to the two quasiparticle disper- 
sions fi ± ^y{W^/2rn^^Jl)^ + A^ and the middle line is that of the 
free atom dispersion. 



D. Momentum Resolved RF 

We turn now to momentum resolved RF spectra. Later 
in Figure [20] we will study some potentially analogous ex- 
periments which address momentum resolved (or "angle re- 
solved") photoemission in the cuprates. In Figure [T3] we 
present an intensity map as a function of single particle en- 
ergy and wave vector k for the (occupied) states which con- 
tribute to the RF current. Here we consider the homogeneous 
case. The yellow to red regions of the figure correspond to 
where the spectral intensity is highest. The temperature here 
is chosen to be relatively high, around 1.9Tc in order to have 
some contribution to the RF cuiTent from pre-existing thermal 
fermionic excitations. The intensity map indicates upward and 
downward dispersing contributions. These coiTespond, to a 
good approximation, to the two RF transitions: to state 3 from 
state 2 with dispersion (i?k + m) and (— i?k + m)' respectively, 
which are shown in Figure [T] These are measured relative to 
the bottom of the band. The width of this contour plot comes 
exclusively from the incoherent terms 7 and So- Here we have 
chosen to represent the latter by an imaginary constant (as in 
the cuprate literature) Sq — —iTq. Based on Eq. ( i26] l and for 
illustrative purposes we take 7 — 0.25Tp and Tq — O.lTp, 
with a small resolution broadening (typical of the experiment) 
as well. 

One can see that the bulk of the current even at this high 
temperature is associated with the pair states which are broken 
in the process of the RF excitation. This figure describes in a 
conceptual way, how this intensity map can be used to com- 
pare with a broadened BCS-like form for the spectral function. 
This form fits very well the two branches shown in the figure 
coiTesponding to upward and downward dispersing curves. In 
this way one can, in principle, establish the presence of pair- 
ing and extract the pairing gap size. 

We stress that these calculations are for the homogenous 
case and it is important to extend them to include the effects 
of a trap. This can be done within an LDA approximation 
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scheme. Once the trap is included the simple analogy between 
the electronic ARPES experiments and momentum resolved 
photoemission spectroscopy is invalidated. However, many of 
the central features survive. While the two branches shown in 
the contour plot in Figure [13] are, in principle present, there is 
a third new branch which appears as well. This corresponds 
to essentially free atoms at the trap edge which will contribute 
significantly |34, 35] to the RF current. It is this branch which 
is also upward dispersing which makes it rather difficult to see 
the effects of the pre-existing thermally broken pairs. 

We summarize the results shown elsewhere ||74ll for the be- 
havior of a unitary trapped gas over a range of different tem- 
peratures. At high T, the central notable feature is a single up- 
ward dispersing curve which fits the free particle dispersion. 
This dispersion can be readily differentiated from that associ- 
ated with pre-existing thermally broken pairs which varies as 
i?k + A* ™d of course, depends on the distribution of energy 
gaps A(r). It arises from free atoms at the trap edge (where 
the gap A(r) is small). As the temperature is decreased to- 
wards Tc a second (downward dispersing) branch becomes 
evident. In the vicinity of the transition, the intensity map is 
bifurcated with two co-existing peaks: one coming from the 
free atoms at the trap edge and the second from the conden- 
sate pairs which are broken in the process of the RF excitation. 
The separation of the two peaks can be difficult to discern un- 
til sufficiently high values of k. Finally, at the lowest tem- 
peratures the striking feature is a single downward dispersing 
branch. This reflects the fact that essentially all atoms are now 
paired in the condensate. Just as in the homogeneous case dis- 
cussed above, a BCS-like fit to this dispersion can be used to 
determine the pairing gap. We stress that there are no abrupt 
changes in the RF behavior at the superfluid transition, very 
much like what we saw earlier in our summary of the cuprate 
literature. 



V. ANALYSIS OF THEORETICAL AND EXPERIMENTAL 
RF SPECTRA 

A. Momentum Resolved Spectroscopy 

We now compare theory and experiment in a trap based on 
the momentum resolved spectra just discussed for the homo- 
geneous case. In Figure [14] we have taken a larger intrinsic 
broadening and included an empirical resolution broadening 
as well, again somewhat larger than the value R w 0.2 indi- 
cated for the experiments. These parameters are seen to opti- 
mize semi-quantitative agreement with the data plotted in the 
top panel from Reference flB] . 

The bottom panel presents the theoretical intensity maps. 
The dotted white curve represents a fit of the experimentally 
deduced peak dispersion while the solid white curve is the 
theoretical counterpart. Here, as in experiment we have fit the 
energy distribution curve to a single Gaussian peak. The com- 
parison between the two white curves shows semi-quantitative 
consistency. Moreover, both the solid and dotted white curves 
can be well fit to the BCS dispersion involving E\^, as was 
originally proposed in Ref. IrTSll . While Figure [T4l seems to 



Contour plot lor Fig. 3b — Spectral function A(k^) l(co) at unitarity and around T=Tc 

1 4500 



t 



I 



3000 
3500 
3000 



\ 




Figure 14: Contour plots of momentum resolved RF spectra in a 
trapped configuration. Top panel (a) is experimental data |20]. The- 
oretical results (b) correspond to occupied spectral intensity map, in 
a unitary trapped Fermi gas at T /Tc = 1.1. Here Eo — Q.i5E% 
and 7 — 0.38-E^ at the trap center. The red curves represent the 
free atom dispersion, while the white solid and dashed curves are the 
quasiparticle dispersion obtained theoretically and experimentally 
1 20], respectively, via fitting the energy distribution curves (EDCs) 
with a single Gaussian. 
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Figure 15: Comparison of calculated RF spectra of a trapped gas 
(solid curve, Tc ~ 0.29) with experiment 13311 (symbols) in a har- 
monic trap calculated at 822 G for two (estimated) temperatures. 
From Ref. (3^ . The dashed lines are a guide to the eye. There is rea- 
sonable agreement, but because final state effects are not included, 
the high frequency tails are overestimated in the theory. 
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Figure 16: RF spectra for a trapped unitary imbalanced gas with 
S = 0.5. Here we take 7 — 0.05. Four different temperatures are 
indicated. Majority spectra are in red and minority in blue. 



Figure 17: Comparison of homogeneous RF spectra of a unitary gas 
with (red,dashed) and without (black) final state effects at T/Tp = 
0.15. The top figure is for 12 superfluid and the bottom is for 13 
case, showing a bound state. Calculations were done with 7 = 6. 



capture the essential results shown in the experiment, with 
higher resolution it should be possible to obtain more direct 
information about the mean experimentally-deduced gap size. 
Importantly, this reasonable agreement and the fact that the 
experiments were done near Tc suggests that there is a sizable 
pseudogap in the Fermi gases at and above Tc at unitarity. 



B. RF Spectra in a Trap 

In Figure [TSl we compare RF spectra in a trap near unitar- 
ity (solid curve) with experiments from Reference ifssll (sym- 
bols) at 822 G on ^Li and for two different temperatures. The 
dashed curve is a fit to the data, serving as a guide to the eye. 
The higher T curve indicates a two peaked structure in both 
data and theory. The lower peak was shown fsst] to be as- 
sociated with nearly free atoms at the edge of the trap. The 
upper peak reflects the existence of pairing. Although it is 
not possible to directly infer the size of the (trap averaged) 
pairing gap A, it is now reasonably clear lfT2|] that a pairing 
gap (pseudogap) is present in the normal state even in these 
early experiments from the Innsbruck group. The lower curve 
can be interpreted to suggest that the atoms at the trap edge 
have temperature small compared to A(r, T). The agreement 
between theory and experiment is not unreasonable for this 
leading order calculation (based on Iq{i')). One can, however 
see that the theory in both cases shows a much slower drop- 
off with increasing high frequency, than seen experimentally. 
We will see shortly that this difference is associated with final 
state effects. 



In Figure [16] we present similar RF spectra for |43] a 
trapped imbalanced gas near unitarity. The polarization is 
S ~ 0.5, and the spectra are plotted for three different tem- 
peratures. Here Tc/Tp = 0.25. It is useful to refer back to 
the lower right panel in Figure |6] to see precisely what region 
of the polarized gas phase diagram is relevant. The red curves 
correspond to the majority and the blue to the minority. For 
the majority, one can see that the free atom peak at — is 
present at all temperatures, unlike the previous case in a bal- 
anced gas. At the highest temperature T/Tp = 0.4 (which 
is close to T* « 0.35) the system is normal and pairing is 
absent. Very close to Tc in the middle panel one sees a clear 
pairing peak signature associated with the pseudogap. For this 
analysis we chose the broadening in Eq. ( |26] |, to be very small 
with 7 = 0.05i?F- 



C. Final State Effects 

As we saw in Figure [15] final state effects are expected to 
cut off the long tails in the RF spectra found theoretically in 
the lowest order theory. It has been argued [24] by the MIT 
group that one should limit the importance of these final state 
contributions by studying a 1 3 superfluid instead near unitar- 
ity with an RF transition involving (for example) 3^2. For 
this case the magnitude of the final state scattering length is 
small, although it is positive and this will lead to a bound state 
contribution in the spectra. One could alternatively argue that 
it is better to work with the 12 superfluid where there are (gen- 
erally) no bound states and where one can, more readily, im- 
pose the sum rules to arrive at estimates of the pairing gap. At 
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Figure 18: Homogeneous RF spectra of unitary, population imbal- 
anced (S = 0.5) gas at very low T = 0.05, where phase separation 
is stable. Upper panel is density profile indicating radii (a) and (b) 
used in lower two panels to compute the tomographic or homoge- 
neous spectra. Red is majority and blue is minority. No final state 
effects are included and we take 7 = 0.05. 



this point both options should be explored. 

In Figure [17] we compare homogeneous spectra at T = 
0.15Tp with and without final state effects for these two dif- 
ferent superfluids at unitarity. The figure on the top corre- 
sponds to the configuration of the Innsbruck ex peri ments [33J 
and on the bottom to recent MIT experiments ll24ll . It can be 
seen that final state effects in both cases do not change iIstIi 
the threshold « 0.4A for the positive continuum (discussed 
in Sectionllllb. However, they do lead to a somewhat sharper 
peak and to a more rapid fall off at high detuning, as is con- 
sistent with the sum rules. 

It is not easy to do the calculations which include final state 
effects in a trap so we can only qualitatively infer from the top 
panel in Figure [IT] that the corrections associated with their 
inclusion are what is needed to improve the agreement be- 
tween theory and experiments on the 12 superfluid in Figure 
[T5] If final states introduced a bound-bound transition on the 
BCS side of resonance, as conjectured [24], it would proba- 
bly not be sufficiently robust with respect to temperature ll74ll . 
and would likely merge with the continuum with increasing T. 
Moreover, this positive continuum contribution, which reflects 
pairing, is always present; theoretically, one never finds sim- 
ply an isolated bound-bound transition. With a small amount 
of lifetime broadening (7) it is likely that a final state induced 
bound-bound transition on the BCS side of resonance would 
not make a noticeable difference. In this way, although some 
concern has been raised |24] about whether or not the Inns- 



bruck experiments were properly interpreted as evidence for 
a pairing gap (rather than a possible bound-bound Ii41i1 transi- 
tion), we concur with their original interpretation. 

By contrast with the 12 superfluid, when final state effects 
for the 1 3 superfluid are included, there is little change in the 
shape of the spectra (shown on the bottom). There is, how- 
ever, an important change of vertical scale associated with the 
bound state. 



D. Tomographic Scans in Imbalanced Gases 

In Figure[T8]we turn to effectively homogeneous spectra as- 
sociated with tomographic plots in the same imbalanced gas 
with S = 0.5 studied earlier Here we consider extremely low 
temperatures so that [48] the system is in the phase separated 
state as can be seen from the lower right plot in Figure |6l One 
can see this phase separation in the top panel which presents 
the density profiles. Also indicated are the two points (a) and 
(b) which establish the radii used in the lower panels at which 
the RF spectra are plotted. Here Rt f — y^2Ep/muj^ is the 
Thomas-Fermi radius and m and ui denote fermion mass and 
trap frequency. The red and blue curves are for the majority 
and minority. The middle panel shows that just inside the su- 
perfluid core there is very little difference in the majority and 
minority spectra as expected for a locally unpolarized super- 
fluid. The lower panel shows the behavior just on the other 
side of the phase separation boundary where there is very lit- 
tle minority and hence very little pairing. As a consequence 
the RF peaks are close to zero detuning. We stress here that 
we are using the strict BCS-Leggett theory without including 
Hartree effects. With the latter included there may be a rela- 
tive shift of the energy scales associated with the majority and 
minority atoms. 

Figure[T9]presents a comparison between theory (with final 
state effects) and experiment | [53l| at moderate temperatures 
(Tc = 0.25Tp) within the Sarma phase. The theory and ex- 
periments are for the 13 superfluid at unitarity. The upper 
panel corresponds to recent data from MIT [53], indicating 
via a contour plot, the various radii probed in the tomographic 
scans. In the lower panel, the counterpart theoretical profile 
indicates the three different radii via (a), (b) and (c). Also 
shown in the theory is where — A would be found within the 
negative continuum. One sees a reasonable correspondence 
between theory and experiment, if we compare b,c and d on 
the top with (a), (b) and (c), respectively on the bottom. We 
stress that Hartree effects have not been included in the the- 
ory so that the zeroes of the horizontal energy scales are not 
equivalent. These Hartree effects have been extensively ana- 
lyzed in Ref. |[53ll , and these authors, as well, have provided 
estimates of the pairing gap. 



VI. PHOTOEMISSION EXPERIMENTS IN THE 
CUPRATES 

We turn now to recent issues in photoemission experiments 
flda in the underdoped cuprates. These call attention to the 
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Figure 19: Comparison between theory and experiment fsS*] in top panel, for tomographic scans of an imbalanced unitary gas. Theory assumes 
5 = 0.5, for definiteness. Here T /Tp = 0.15. The calculations are for a 13 superfluid (as in experiment) with final state effects included 
(1 /kpau = 2.5), and 7 = 0. Red indicates majority and blue, minority. Panels a, b and c in the theory should be compared with b, c and d in 
experiment. 



question of how the k dependence of the spectral function 
varies as one crosses Tc- The earlier discussion in Figure [3] 
focused on the k integrated photoemission spectra. As in this 
previous figure here we address how superconducting coher- 
ence is manifested when there is a normal state pseudogap. 
We stress that measurements like photoemission and RF spec- 
troscopy are not phase sensitive probes of the system and can- 
not directly prove the existence of superfluidity. 

Here we focus on the regime near the gap nodes, where the 
gap is smallest. The experiments of interest here very likely 
contain important clues as to the nature of the superconduct- 
ing state which appears in the presence of a normal state pseu- 
dogap. While many aspects of the cuprates below Tc appear 
to be typical of (c?-wave) BCS superconductors one expects 
some differences to appear simply because an excitation gap 
is present at the onset of superconductivity. 

There are four key points which have been identified in 
Ref. It^. We believe these are consistent with a BCS-BEC 
crossover interpretation of these cuprate photoemission data. 
As reported IITql (i) the excitation gap, A(k), as measured in 
photoemission experiments remains roughly constant in tem- 
perature from very low T, to temperatures well above Tc- 
(ii) In the superfluid phase A(k) displays the expected point 
nodes (associated with d-wave symmetry); however, these 
rapidly broaden into Fermi arcs once the temperature reaches 
the vicinity of, and surpasses Tc- Importantly, "this remark- 
able change occurs within the width of the resistive transition 
at Tc". (iii) It has also been reported iIttIi that the energy scale 
associated with the excitation gap appears to be T*, which 
is conventionally taken as the pseudogap f ill [T2I [TtIi onset 
temperature, and that the Fermi arc length scales with T/T* 
above Tc- From (i) it is inferred that (iv) "the energy gap is 
not directly related to the superconducting order parameter". 



To address these and other photoemission experiments, the 
normal state self energy is taken ||63[ l64ll66|] to be of the form 
shown in Eq. ( |26] | with Eq. ( l23T l and with a purely imaginary 
background self energy: I]o(k, = — iFo. The rapid, but 
smooth destruction of the d-wave point nodes can be phys- 
ically associated with the fact that the superconducting or- 
der parameter A^c disappears smoothly but precisely at r+. 
Above Tc the effects of 7 and Eq lead to a smearing and the 
point nodes are replaced by Fermi arcs |63]. Below Tc with 
the onset of phase coherence through A^c, the arcs are rapidly 
replaced by point nodes. One says that there has been "a col- 
lapse of the Fermi arcs", and that the nodes are "protected" 
below Tc. 

The collapse is a continuous process. We argue that it is not 
to be associated with an abrupt disappearance of the inverse 
lifetime 7, but rather it reflects the gradual emergence of the 
condensate to which the finite momentum pairs are continu- 
ously converted as T decreases. This is related to the fact that, 
from Eq. ( |23] ). we see there are two terms in the self energy 
below Tc- At the lowest temperatures Apg vanishes, whereas 
above Tc, it follows that Age is zero. 

In Figure|20]we address these new experiments by showing 
the collapse of the Fermi arcs from above to below Tc within 
our BCS-BEC crossover formalism; we plot the percentage 
of Fermi arc length as a function of T/T* and for different 
doping concentrations from the optimal (T* k, Tc) to the un- 
derdoped regime (T* >> Tc). The observed collapse is inti- 
mately connected with our earlier observation that the spectral 
function in Eq. ( |27] i has a zero ?Au} — — below Tc, whereas 
the spectral function has no zero above Tc- 

In summary, this reasonable agreement between theory and 
experiment shown in the plot supports our physical picture 
that pseudogap effects persist below Tc in the form of non- 
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Figure 20: (color online) Fermi arc length as a function of T/Tex for 
doping concentrations from optimal to underdoping for a cuprate su- 
perconductor. Fermi arc length is typically finite above Tc and drops 
to zero upon the onset of phase coherence. The normal state portions 
of the curves is close to universal, in agreement with Ref. [TJ. The 
comparison in the inset between the theory with a slightly (15%) 
enlarged Tex and experimental data (symbols) |76] shows a good 
semi-quantitative agreement. 



condensed pair excitations of the condensate. We argue that 
the pre-formed pairs above Tc do not abruptly disappear just 
below Tc- This recognition that the superfluid phase is, it- 
self, very complex represents a new direction for the cuprate 
field which these important ARPES experiments have now 
called to our attention. We also note that this purportedly 
more complex superfluid, containing non-condensed pair ex- 
citations can potentially be elucidated through studies of BCS- 
BEC crossover in the cold Fermi gases. 



VII. CONCLUSIONS 

There has been enormous progress in the field of radio fre- 
quency (RF) spectroscopy of the Fermi gases. This technique 
holds the promise of being as valuable to these atomic super- 
fluids as photoemission has been to the cuprates. We have 
tried in this Review to argue that it also holds the promise of 
helping to address some fundamental issues in the cuprates 
which are very general, such as how to describe that anoma- 
lous superfluid phase which forms in the presence of a normal 
state excitation gap. 

On a less general level, these RF experiments also hold the 
promise of helping to address (that is, support or rule out) 
one particular approach to the theory of high temperature su- 
perconductivity: namely that based on BCS-BEC crossover. 
There are many alternative physical pictures of the cuprates 
and, indeed, the ultracold gases on optical lattices have pre- 
sented themselves as possible simulators of at least some of 
these alternatives. Most notable among these are proposed 
cold gas studies of the repulsive Hubbard model which is 
thought to be relevant to the "Mott physics" aspects of the 
high Tc superconductors. 

A goal of this Review was to present a broad background 



to the theory of RF spectroscopy and its relation to photoe- 
mission spectroscopy. We summarize key experiments us- 
ing both techniques and show how they can be addressed 
within a BCS-BEC crossover approach. Included in our anal- 
ysis are trapped as well as homogeneous gases, and popula- 
tion imbalanced gases. The RF field has seen a proliferation 
of methodologies including tomographic and momentum re- 
solved scans, all of which are discussed here. 



The immediate excitement surrounding these RF experi- 
ments is in large part because they hold the potential for mea- 
suring the pairing gap A. As time passes, however, it has be- 
come progressively more clear that extracting detailed quan- 
titative information is increasingly difficult. There is no clear 
signature in the RF spectra of A except to establish whether 
it is present or not. This observation applies to both balanced 
and imbalanced gases. Here we summarize two methodolo- 
gies which might make it possible to extract more precise 
numerical values for A: either through sum rules which in- 
volve integrals of the measured RF spectra, or through recent 
momentum resolved experiments. At this moment, however, 
neither of these has been put into experimental practice at a 
quantitative level. 



Our analysis (and all results presented in this paper) are 
subject to the caveat that they are obtained within the BCS- 
Leggett crossover theory (extended to finite T). This theory, 
like all others is not an exact theory. However, this approach is 
one which is particularly well suited to the cuprates because 
there are no spurious first order transitions. Recently, there 
have been RF calculations based on an alternative scheme 
which is closer to the Nozieres Schmitt-Rink approach to 
BCS-BEC crossover. This is very nice work, which we only 
have only briefly mentioned, which addresses radio frequency 
experiments in the normal state and at T w 0. The latter, in 
particular has been important in elucidating how to incorpo- 
rate final state effects. We also only make passing reference 
to the body of work on the highly imbalanced gases. 

It should be clear that this field is moving rapidly and much 
has been accomplished at the theoretical and experimental 
ends. However, in writing this Review, we felt it was timely 
to suggest possibly new directions for cold gas research which 
might elucidate the cuprate superconductors. As summarized 
earlier, the key issues which have emerged in photoemission 
studies involve characterization of the fermionic self energy, 
of the pseudogap and of the signatures of superconducting co- 
herence (in passing from above to below Tc). These issues 
have a counterpart in the cold Fermi gases as we have sug- 
gested here and it will be of enormous benefit in future to 
exploit this new class of "materials" to address these funda- 
mental questions in condensed matter. 
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Figure 21: RF current of the unitary polarized {5 — 0.1) gas at T = 
O.ISTf and T = 0.25rF. For these temperatures, the Sarma state is 
stable. The black curve is for the majority and red is for the minority. 
The arrows indicate the position of — A. Tc ~ 0.24rF in this case. 
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Figure 22: RF current of the unitary unpolarized gas at T = 0.15Tf 
and T = 0.25Tf . = 0.255rF in this case. This figure is included 
as a contrast with the polarized case, to make it evident that the nega- 
tive continuum is more difficult to see in the absence of polarization. 



Appendix A: ANALYTICAL RESULTS FOR RF SPECTRA IN 
HOMOGENEOUS, POLARIZED GAS 

Following the same arguments as in Section III El for the 
polarized case, we have for the majority current 



^^ko[l- f{Eo + h)], (i.>i/i)(Al) 

= ^^kof{Eo-h), < < 0) (A2) 
and minority RF current 

^o'^(^) = ^^^o[l-f{Eo-h)], (z.>:/i)(A3) 



= ^^''of{Eo + h), < < 0) (A4) 

Here /i — + /i|)/2 is the average chemical potential, and 
/i = /i| — /i| is chemical potential difference. Since Eq + h 
is always positive, the minority RF current is not associated 
with a negative continuum at low temperatures. On the other 
hand, since generally Eq — h < 0, the majority RF current 
has a negative continuum. For strictly zero temperature, the 
negative cont inuum is located in the range —h — y/h'^ ~ < 
V < —h + ^/h?~~A? . In Figure |2T1 we show the RF cuiTent 
for the stable Sarma state for small population imbalance 5 = 
0.1. We note that the negative peak location is not at — A. 

In Figure|22]we compare with the unpolarized case, show- 
ing that there is essentially no negative continuum (until one 
approaches temperatures closer to Tc). We can ask why the 
onset of the negative continuum is unrelated to the energy 
scale A. The answer is that while, the negative continuum 
in the Sarma phase is associated with the presence of excess 
fermions (as seen by comparing the two figures in this ap- 
pendix), these excess fermions in the Sarma phase are essen- 
tially gapless. Thus, they are not directly associated with an 
energy scale which involves A in any explicit way. 
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Appendix B: DETAILS ON THE FINAL STATE EFFECT 
DIAGRAMS 



The response function which yields the RF current is 



D{y)=Do{v) + Dl[y)t,^{y) 



(Bl) 



where the first term is the usual leading order contribution and 
the second corresponds to the Aslamazov-Larkin (AL) dia- 
gram given by 



Dal[Q) = - G^^HK)Gi'\-K)G^^\K + Q) 

K.K'.P 



AY.G{K)Go{-K)Gf{K + Q) 



K 



where, as discussed in the text we have ~ — A^/T. This 
term can be rewritten in terms of the functions F{Q), and 
D2Q defined in the main text. We introduce the RF detuning 
frequency ly 



(B2) 



The retarded form of the various complex functions which 
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enter into the RF current can be summarized as 



/(i?k) - /(Ck,3) 2 



1 - /(i?k) - /(a^s) 2 



- £^k - Ck + ^0^ 
A rl-/(£;k)-/(ek,3) 



2-Bk 



- £^k - Ck 



^0^ 



/(i^k) - /(^k,3) 



^0^ 



1 - /(i?k) - /(^k.s) 



/(i^k) - /(^k. 



^0^ 



2ek 



(B3) 



(B4) 



(B5) 



In the general case with ai2 7^ 013, it follows that 

ImD^(z^) = ImD(f 

-Re(J'^)^Im^^3^--" + ImjD^)'^ Re^rs^"" 



After a straightforward calculation we have 



A / TO TO 

Numerator ~ — C 



V 47rai3 47rai2 



(B14) 



In this way we rewrite the response function in the form 
presented in the text as 



f m m \ lmDQ{i') 



V47rai3 47rai2/ \t^rl'^{v 



(B15) 



To compute the RF current, we calculate the real and imag- 
inary parts of the above quantities. For later convenience, we 
introduce the following expressions. 



E — 



1 - fiE^,) - /(a,3) /(i?k) - /(a,3) 



- -Ek - Ck + -Ek - ^k 

K 

(B6) 

C{u) = _7rE^{[l-/(^k)-/(^k,3Mj^-i?k"a) 



Appendix C: SPECIAL CASE OF EQUAL INTERACTIONS 



In the special case ai2 = 013, we expect I{u) oc <5(i/). 
We showed above that Imtj^j^' = —Cv — _B 0+. The second 
term plays essentially no role when ai2 7^ 013 case, but it is 
important in this special case. Combined with the real part, 
we have 



tl Q 



[/(i?k)-/(ek.3)]'5(i' + i?k-a)} 



(B7) 



We also define = 0) ee B{0) and C{v = 0) = C(G). For 
the most part it will be convenient to drop the argument and 
simply write B, C. 
It can be shown that 



ImDi 



ReD^ 



A2 
h' 

AC 
AB 



C 



(B8) 

(B9) 
(BIO) 



We notice that 



'E — 

^ 2£;k 



1 - /(i?k) - /(a,3) f{Ek) - fi^Ks) 



- -Ek - Ck 



{B + iC)i^ + iBO+ 
Using the fact that C — (when ai2 = ai3) we find 

IB 1 C TT „ , 

^-^WTc^' '''''' = -.WTc^ + B^^"^ 

Thus we have 

Imi:>^(i^) = ImL>^(iy) +Re£l2lmtf3 +ImD2""Retf3 
= ttA^BS{v). (CI) 

In the limit 012 = 013 

A2i?(o) = -E{[i-/(^k)-/(ek,3K 

k 

+ [/(i^k) - /(a.3)]^} 

= -in2 - n^) . (C2) 



-Bv 



i' + Ek-(,k 

(BID 



Thus we deduce 



1 



/(i.) = --lmD"{u) = (712 - n3)S{iy) , (C3) 



We make use of the 1 , 2 superfluid gap equation and identify 

Xi2(0) - Xi3(0), so that 



-1,R 



TO 



13 



47rai3 47rai2 

TO TO 



Imt 



-1,R 



13 



47rai3 47rai2 

Im[xf3H-Xi3(0)] 
-Cv-BO+ 



+ X13M - Xl3(0) 
- Bi' (B12) 



(B13) 



which is the expected II37I1 result. 



Appendix D: DETAILS ON EXPLICIT EVALUATION OF 
SUM RULES 



We will next prove the sum rules on the zeroth and first 
moment of the RF spectra for the unpolarized case. We focus 
on the general case 012 ^ 013. In the large limit, from 
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Eqs. dlSll and dBT) it follows that B - i^'^/^ and C - z^'^/^ 
Thus we have Xisi'^) ~ Xi3(0) ^ Biy ~ v^^^ as — > oo. We 
deduce ti3(i^) = / gis + Xisi^)]^^ v^^l'^. Thus we have 
/o ~ v^^l"^ and / ~ i/^^/^ 2&v^oo. Consequently, we see 
that final state effects convert a diverging first moment into a 
finite result for the clock shift. 

To obtain this shift quantitatively, we use the Kramers- 
Kronig relations for iia. 



For convenience we define (5g — 



TT V — V' 



(Dl) 



47rai3 4x012 

bg-Bv 



so that 



(bg - BvY + (Cz/)2 



Then from the Kramers-Kronig relations with = we 
have 

We also have by taking a derivative and setting = 0, 
From these results we can prove the sum rule 



2.2 / dv\mt^j^{v) 



TT V 

-B(0)A' = ^2-^3 



(D4) 



Then the first moment is 



y dvvliy) = -Sg^A^ J ^ 



2 2/" dvlmt^^iiy) 



m / 1 1 



= -5gA' = AV 

An \ai2 013/ 

Combining these equations we find for the mean clock shift 

A2 m / 1 1 



J dh' I{h') 712 — ^3 47r \ai2 ai3 



(D6) 



which is the expected result [37]. In the text we explore how 
this sum rule may be used to measure the pairing gap. 



Appendix E: FINAL STATE EFFECTS IN A 
HOMOGENEOUS BUT POLARIZED SYSTEM 



In the polarized case, state 1 and 2 have different chem- 
ical potentials. All the calculations will closely parallel the 
unpolarized case. In the following, we only consider the RF 
transition from state 2 to 3. The formulae for the transition 



from from state 1 to 3 can be similarly derived. As before, we 
define 

F'{K) = -AG'^^\K)Gi^\-K) 

1 —iuji + + h 



-A- 



ii^i - h {iLUi - h) 
A 



--El 



[ioji - hY - El 



(El) 



Here pi = (^^ + ^i^)/2 and ft = (^^ - /iJ/2, 
Generalizing D2{Q) to the polarized case we find 



k " 



1 - /(gk.T) - /(ek.3) 
(irj„ + h) - Ei^- ^k,3 

) - /(ek,3) 



{i^n + h) + - ek.3 



with i?k,T — E\^ — h and i?k,x — E\f^ + h. 

We analytically continue along with a shift of variables to 
write: 



(E2) 



In this way we find the denominator in D2{Q) is the same as 
for the unpolarized case; the only difference is in the numera- 
tor. 



1 - /(i?k) - /(ek,3) 
/(i?k) - /(Ck,3) 



l-/(i?k,T)-/(^k,3) 
/(i?k,i)-/(ek,3) 



We can similarly evaluate Do{v) and Xi^i^)- Finally, we 
not that 



with 



Bo 





A2 

—C2 

v 






(E3) 




m 


m 


B2iy 


(E4) 


47rai3 


47rai2 


Imi 13 


-C2iy 






(E5) 


^ 1 ri- 


f{Ew.l) 


- /(fk,3) 


fiEKl 


- /(ek,3) 



2E^ 



I' + Sk - a 



C2M = -^^^{[i-/(i?k,T)-/(Ck,3M^^-i?k-a) 

k 

- [/(i?k,i)~/(a,3)]% + £^k-a)} 

In this way we find for the RF current of state 2. 
1 / m 



which is formally the same result we found in the absence of 
population imbalance. These results can readily be general- 
ized to compute the current in state 1 . 
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1. Sum Rules for the Polarized Case 



RF current 



Just as for the unpolarized case, we have the Kramers- 
Kronig relations (IDll i for tf^. Following the same analysis 
as for the unpolarized case we have 



+ 00 



l/Sg 



Now using the fact that A^_B2(0) — —{n2 — ns), it is 
straightforward to obtain the zeroth and first moment of the 



Then the average clock shift is 



/(it/ 12(1^) n2-n3ATr\ai2 013^ 
Again, this is formally, the same as for the unpolarized case. 
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